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ABSTRACT 

We present estimates of the gravitational lensing shear variance obtained from images taken at the CFHT using 
the UH8K CCD mosaic camera. Six fields were observed for a total of 1 hour each in V and I, resulting in catalogs 
containing ~ 20,000 galaxies per field, with properly calibrated and optimally weighted shear estimates. These 
were averaged in cells of sizes ranging from 1'.875 to 30' to obtain estimates of the cosmic shear variance (7 2 ), 
with uncertainty estimated from the scatter among the estimates for the 6 fields. Our most reliable estimator for 
cosmic shear is provided by the cross-correlat ion of the shear measured in the two passbands. At scales 10' 
the results are in good agreement with those of Van Waerbeke et al. (200C ), Bacon, Refregier, & Ellis (2000 ) and 



Wittman et al. (2000| ) and with currently fashionable cosmological models. At larger scales the shear variance 



falls below the theoretical predictions, and on the largest scales we find a null detection of shear variance averaged 
in 30' cells of (j 2 ) = (0.28 ± 1 .84) x 10" 5 . 

Subject headings: Cosmology: observations — dark matter — gravitational lensing — large-scale structure of 
Universe — galaxies: photometry 



1. INTRODUCTION 



Weak lensing provides a potentially powerful probe of mass 



fluctuations in the Universe (Gunn 1967; Mellier 1999 and ref- 



erences therein). Three independent groups have recently pre- 
sented estimates of the s hear variance from deep 'b lank-field' 
CCD imaging surveys. Van Waerbeke et al. 2000] (hereafter 
vWME+) measured the s hear variance in circular cells o f radii 
ranging from 0'.7 to 3'. 5; Bacon, Refregier, & Ellis 2000 (here- 



after BR E) measured the shea r variance in square cells of side 
8'.0 and [Wittman et al. 2000| (hereafter WTK+) have provided 
estimates of the shear-shear correlation function at separations 
3'. 25, 8'. 5 and 22'. 0. Here we present shear variance measure- 
ments from ~ 1.5 square degrees of deep photometry obtained 
as part of our ongoing weak lensing survey. We find results 
which are broadly in good agreement with the recently pub- 
lished estimates. 

2. the data 

The data were taken at the 3.6m CFHT telescope using the 
8192 x 8192 pixel UH8K camera at prime focus. The camera 
delivers a field size of 0°.5 with 0".207 pixels. Our survey strat- 
egy has been to target blank fields in six widely separated areas 
for ease of scheduling, and in each area we plan to make 6 or 
so pointings scattered over a region of extent ^3°. In January 
1999 the UH8K was replaced by the CFH12K camera. By that 
time we had completed 1 hour integrations in both V and I for 
two pointings in each of three areas. The field names, centers 
and also the estimated seeing are given in table [l]. One of the 8 
devices in the mosaic (lying in the NW corner of our images) 
has very poor charge transfer efficiency and the data from this 
device were discarded. After further masking of regions around 
bright stars the total useful solid angle per field is ~ 0. 16 square 
degrees. 

3. DATA REDUCTION 

The data w ere reduced much as for our MS0302 supercluster 
observations (Kaiser et al. 1999). After flat fielding, an ob- 



ject finder was applied to each image, and a set of bright but 
non-saturated stars extracted for registration purposes. The po- 
sitions of these stars, a long with celestial coordinates from the 
USNOA (Monet 1998) catalog, were used to find a mapping 



from image pixel coordinates to orthographic sky coordinates. 

Images of the stars were analyzed to generate a model 
for the point spread function (PSF) g(x;r), this being the 2- 
dimensional profile of a star with centroid at r measured in 
coordinates x being measured relative to the centroid. The 
model is a sum of 2-D image valued modes g,(x): g(x;r) = 
J2i c i( r )gi( x ) with coefficients c,(r) which are low order poly- 
nomials in star position. We found a 1st order model to be 
adequate to describe the variation of the PSF with position on 
the chip (though see below). 

The astrometric solutions, the PSF models, and also standard 
star observations were first used to create a set of photometri- 
cally calibrated 'homogenized' images which were degraded to 
have a common identical PSF. These images were compared in 
order to identify cosmic rays and other transient events. Next, 
for each raw image we generated a re-circularized image by 
convolving with a 90-degree rotated version of the PSF model. 
These, as well as the raw images, were then warped to sky coor- 
dinates, with previously identified cosmic rays being removed, 
and the stacks of images combined to provide a quilt of overlap- 
ping images. This procedure results in three images: a median 
of the raw images (in which the PSF is generally non-circular), 
a median of the re-circularized images, and also a sky-noise im- 
age. The final summed images were sampled with 0".15 pixel 
size. 

The final object catalogs were obtained by applying 
hf indpeaks to the median averaged raw images, this pro- 
gram having been modified in order to allow properly for the 
non-trivial noise correlations in the final images. After apply- 
ing aperture photometry analysis, shape s of the objects were 
measured as described in Kaiser (2000| ). The essential result 
of this is a polarization vector q a = M a i,„ J d 2 x w(x)x;x,„/(x) 
formed as a combination of weighted second moments of the re- 
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Table 1 
Field centers and seeing 



area 


pointing 


RA (J2000) 


DEC (J2000) 


/ 


b 


FWHM(I) 


FWHM(V) 


1650 


1 


16:51:49.0 


34:55:2.0 


57.37 


38.67 


0".82 


0".85 




3 


16:56:0.0 


35:45:0.0 


58.58 


37.95 


0".85 


0".72 


Groth 


1 


14:16:46.0 


52:30:12.0 


96.60 


60.04 


0".80 


0".93 




3 


14:9:0.0 


51:30:0.0 


97.19 


61.57 


0".70 


0".85 


Lockman 


1 


10:52:43.0 


57:28:48.0 


149.28 


53.15 


0".83 


0".85 




2 


10:56:43.0 


58:28:48.0 


147.47 


52.83 


0".84 


0".86 



circularized image, and a polarizability tensor which describes 
the response of the polarization to gravitational shear. The 
weight function w(x) was taken to be a Gaussian ball of 2 pixels 
in scale length. For convenience, the quantities actually gener- 
ated were the normalized polarization q a = q a j y/qpqp and a 
polarizability Q defined such that for galaxies with this shape 
and size, the expectation value of q a is (q a ) = Qja- The q a val- 
ues were corrected for artificial shear introduced in the image 
warping by slight errors in our astrometric solution. 

Finally a selection on significance of 4 < v < 100 was ap- 
plied to select a 'faint galaxy' catalog. The corresponding mag- 
nitude limits are somewhat fuzzy since significance depends on 
both size and magnitude. The counts (number per magnitude 
interval) turn over at about ni\ ~ 24 and my ~ 25. Numbers of 
objects in the final catalogs are shown in table ||. The density of 
objects on the sky is very similar to that in the images obtained 
by vWME+, WTK+, and slightly higher than the density of ob- 
jects in the BRE sample, so our shear variance estimates should 
be more or less directly comparable. 

Given a single galaxy, a fair (but very noisy) estimate of the 
shear is j a = q a /Q- However, since the normalized polariza- 
tion response Q varies with shape and size of the object (small 
objects having very little response for example), to measure the 
mean shear in a region containing N galaxies one should weight 
the individual estimates by Q 2 so the optimal mean shear esti- 
mate is 7 Q = Xi2?a/E2 2 - This assumes that one has little 
prior knowledge of the of galaxy redshifts. If the averaging re- 
gion contains a large number of galaxies, as is the case for the 
cells considered here, one can replace Q 2 by N(Q 2 ) where 
the (Q 2 ) is an average over all of the galaxies in the catalog. 
The optimal mean shear is then the average of weighted shear 
values for the individual galaxies: jc = ^cqc/Qc, with nor- 
malized weight ujQ = Qq/(Q 2 )- The quantity (Q 2 ) is a useful 
measure of the image quality. It is equal to the inverse shear 
variance per galaxy, so, for instance, the statistical uncertainty 
in the mean shear 7 measured from a sample of N galaxies is 
a 2 . = 1 /(N(Q 2 )). The values of (Q 2 ) are also given in table |. 

Prelimi nary results of this analysis (W ilson, Kaiser, & Lup- 
pino 1999) in the form of estimates of the net shear for each 
of the six fields and for catalogs generated from the I and V 
images separately gave shear values typically of about 1%, but 
with a few larger values. These, however, showed little corre- 
lation between the two passbands, suggesting that the results 
were contaminated by some systematic error. Examination of 
the re-circularized images of the stars in the fields with seem- 
ingly spurious shear values revealed at least a major part of the 
problem. The mean stellar polarization was found to vary sys- 
tematically with magnitude. This is to be expected for very 
bright stars where the pixels saturate and charge begins to bleed 



along the slow direction of the CCD. The effect found here had 
the same signature (a trend for q\ to become negative for bright 
objects) but appeared at a low level and, unexpectedly, for stars 
much fainter than the saturation limit. The result was that in 
some cases our PSF model fitting procedure, which weighted 
stars according to brightness, did not correctly recircularize the 
faint stars as it should (since we require the PSF appropriate 
for faint galaxies). The effect seemed to be variable, and also 
tended to be associated with particular chips. As a simple fix, 
we fit the residual polarizations for the faintest stars (/ > 18, 
V > 21) to a 4th order spatial polynomial and then used a smear 
polarizability ana logous to that defined by K aiser, Squires, & 
Broadhurst (1995) to correct the galaxy q a values. This re- 
duced the spurious shear values considerably. 

Table 2 
Galaxy catalogs 



field 


Ni 


<e 2 )i 


AV 


(Q z )y 


1650 1 


21569 


1.527 


15403 


1.234 


1650 3 


18187 


1.785 


16518 


1.147 


Groth 1 


27293 


1.440 


16391 


1.543 


Groth 3 


19162 


1.254 


15876 


1.775 


Lockman 1 


20726 


1.855 


20358 


1.352 


Lockman 2 


20017 


1.630 


17779 


1.417 



4. COSMIC SHEAR VARIANCE 

We have chosen to focus here on a single simple statistic: the 
variance of the shear averaged in cells of various sizes. This is 
the statistic used by vWME+ and BRE and is simply related to 
the shear covariance function presented by WTK+. The cell 
averaged shear variance is also si mply computa ble from the 
spectrum of mass fluctuations (e.g. Kaiser 1992 ), so this pro- 
vides a useful link between observation and theory. Now each 
weighted shear estimate consists of a random intrinsic com- 
ponent 7ci nt and a 'cosmic' component proportional to the in- 
tegral of the tidal field along the line of sight. Modeling the 
cosmic shear as the sum over a set of statistically independent 
screens, we have 



7C = 7Gint + U>G*^2 ls{dc)P(ZG , Zs) 
S 



(1) 



where 7s(0) is the shear field for the 5th screen and for ficti- 
tious sources at infinite distance, and (3 = max(0, 1 -Dsc/Doc) 
is the usual ratio of angular diameter distances. 

This model allows one to compute the variance of the shear 
averaged over galaxies falling in a cell on the sky. We will also 
be interested in the co-variance of shear measured in different 
passbands. Consider the mean shear j P = (l/Np)J2lP f° r a 
specific cell and for galaxies found in two passbands P =A,B. 
Averaging over realizations of random intrinsic shear values, 
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and also averaging over an ensemble of realizations of cosmic 
shear screens, yields 

fa -^ = WW fa • ^ + EW«) C$ 

where iV^s is the number of objects in the cell which were de- 
tected in both passbands. This formula is also valid when A and 
B are the same. The expectation value of the dot product of the 
cell averaged shear is therefore equal to a noise term plus a cos- 
mic term which is a sum of the cell-averaged shear variances 
for the screens. Interestingly, the noise term involves the total 
shear variance (*y A -j B ), containing both intrinsic and cosmic 
contributions, which is convenient since this is the quantity that 
one can actually measure. In obtaining (Q) we assumed that 
the faint galaxies are randomly distributed on the sky, this be- 
ing motivated by the fact that the angular correlation function 
is very small, with w(9) <, 10~ 2 on all relevant scales. 

To implement this, for each field we averaged the shear in a 
grid of contiguous square cells of side L, and those cells in the 
lower quartile of occupation number were discarded. To obtain 
an estimate of the cosmic shear variance we then computed for 
each field 



(7 2 ) ab = — fa ■ 7b- Nab (l A ' 1b) /(N a N b )] 



"cells 



(3) 



cells 



The shear covariance functions are • 7j) = l/(2 2 )i an d 
(Tv ' 7v) = 1/ (G 2 )v whereas (7, ■ 7 V ) was estimated by cor- 
relating the shears for objects which were detected in both the 
I and V catalogs. The shear variances estimated from the sepa- 
rate fields were then averaged together to obtain a final cosmic 
shear variance with uncertainty estimated from the scatter of 
the field estimates about the mean. 

The diagonal components of (t" 2 )ab provide estimates for 
the shear variance for the respective passbands, with strength 
roughly proportional to the square of the mean distance to the 
galaxies (assuming a spectrum of mass fluctuations with in- 
dex n ~ — 1), and the off-diagonal components should lie some- 
where in between. The results are shown in figure |]and in table 
|3| and deviate somewhat from this expectation: the I-V cross 
correlation lies systematically below both the I- and V-band 
shear variance estimates. These results are robust to changes 
in the order of the polynomial in the stellar polarization model, 
and they are not caused by a few discrepant cells. 

The difference between the I and V band shear variance 
may be due to differences in the r edshift distributions, some 
evidence for which was found by Luppino & Kaiser (1997) 
in their study of MS 1054. However, we typically find about 
60% of the galaxies are detected in both passbands, so this re- 
quires fairly high redshifts for the blue galaxies. For example, 
assume that the I-band sample has a redshift distribution like 
that measured by Cowie (personal communication) in the range 
23 < Iab < 24.0, but that the V band sample contains an addi- 
tional 40% population of higher redshift galaxies. If we place 
these at redshift 3, we find that the V-band shear variance is 
about a factor 2 higher than the I-band, much as seen in figure 
[|. However the I-V cross-correlation is then predicted to be 
about 30% higher than the I-band variance, which is not seen. 

The simplest interpretation of these results is that the shear 
inferred from the I- and V-band data separately has been in- 
flated by residual systematic errors of some kind. The level of 
these errors is on the order of 2 percent rms shear on scales of a 



few arc-minutes, falling to somewhat below the 1 percent level 
on 30' scales. If so, the most reliable estimate of the cosmic 
shear variance is provided by the I-V cross-correlation since 
systematic errors which are uncorrelated between the passbands 
will cancel out. Of course there is no guarantee that the cross- 
correlation is not affected by some source of error which is 
common to both passbands, a specific example of which is ar- 
tificial shear arising from intrinsic correlation of galaxy shapes 
in clusters etc. due to tidal effects. 
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FIG. 1 . — Estimates of the cosmic shear variance as a function of averag- 
ing scale L. Lower panel shows the same data multiplied by L to show more 
clearly the large angular scale results. The dotted lines are predictions for a 
currently fashionable cosrookigica] mnHpl — duster normalised; f2 m = 0.3; 
Ca = 0.7; r = 0.25 — from [Tain & Seljak (1997| ) for sources with an effective 
redshift (i.e. equivalent single screen redshift) of z e ff =1-0 and z c ff = 2.0. The 
points have been displaced slightly laterally for clarity. 



Table 3 
Shear variance 



L/arcmin 


10 4 x ) n 


10 4 x (^>vv 


10 4 x (7 z )iv 


30 


0.15±0.18 


0.59 ±0.32 


0.03±0.18 


15 


0.87 ±0.42 


1.57±0.54 


0.48 ±0.22 


7.5 


2.06 ±0.75 


4.77 ±0.83 


1.48 ±0.93 


3.75 


3.51 ±1.40 


9.69 ±2.09 


2.44 ±0.92 


1.875 


5.66 ± 1.40 


11.38 ±2,29 


4.06 ± 1.57 
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FIG. 2. — Our estimates of the cosmic shear variance from the I-V cross- 
correlator are shown as the heavy points. Also shown are results from 
vWME+, BRE, WTK+. The error bars on the vWME+ estimates are statisti- 
cal only. All others are total error including cosmic variance. The lower panel 
again shows the same data multiplied by L t o, show more clearly the large an- 
gular scale results. The dotted lines are the Iain & Seljak (1997) predictions 
as in figure hi 

The I-V shear variance estimator is shown with an expanded 
vertical scale in figure |[ Also shown are the recently an- 
nounced results. The BRE result is shown as presented in 
their paper and with total error estimate including cosmic vari- 
ance. The vWME+ circular cell average shear are plotted 
against L = y/irO. The vWME+ error bars are statistical only. 
WTK+ presented estimates of the ellipticity correlation func- 



tion Ci(ff) = (ei(O)ei(f?)). We have converted th eir C\(9) to an 
equivalent shear variance using formulae from Kaiser (1992) 
with e = 27 and assuming a spectral index n = — 1. The lower 
panel shows the variance multiplied by averaging box size L. 
For a n = -1 spectrum, corresponding to a mass auto-correlation 
function oc 1/r 2 , this quantity should independent of scale. 

At small scales < 10 arcmin there seems to be remarkably 
good agreement between the independent estimates. Note that 
the measurements were made using three separate observing 
facilities. At L = 3'.75 we find ~ 2.5 ~ 10" 4 . This about 
a factor 4-5 lower than the prediction for a light-traces mass 
tt m = 1 cosmolo gy, and an eff e ctive redshift for the background 
galaxies z e ff = 1 ( [Kaiser 1992| ; [lain & Seljakl997| ). 

At larger scales the shear variance we find falls below that 
of WTK+. Their largest scale estimates appear to conflict with 
our null result at about the 2-sigma level. Our large-angle re- 
sults are also smaller than the f2 m = 0.3, ft\ = 0.7 theoretical 
model predictions. 

5. DISCUSSION 

For an effective background galaxy redshift of z e s —1.0 
these measurements probe mass fluctuations in a shell peaked 
at z — 0.4. At this redshift the 30' field size corresponds to a 
comoving distance of about 6/r'Mpc, so the cell variances pre- 
sented here probe scales in the range 0.4-6/z _1 Mpc. On the 
smaller end of this scale we find very good agreement with re- 
cently announced estimates from other groups, and also with 
canonical cosmological theory predictions. It is hard to defini- 
tively rule out the possibility that the small angle measurements 
are inflated by systematic errors, but one can safely rule out 
theories such as light-traces mass high density models which 
predict shear variance a factor ~ 5 higher than our results. 

On larger scales our measurements are extremely precise, 
yet we find only a null detection for our largest cells. These 
results show that on large scales the rms shear is at most a 
fraction of a percent. The apparent discrepancy between these 
results and the theoretical predictions is quite interesting, and 
suggests a steepening of the mass correlation function at scales 
~ 1 -2/r'Mpc. More data are needed however to definitively 
confirm this. 
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